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INTRODUCTION 


Proper  biomechanical  descriptions  of  major  human  joints  are  essential 
for  improving  the  biodynamic  response  of  the  multisegmented  models  of  the 
human  body.  In  a  series  of  technical  reports  [1-3]  and  articles  [4-18] 
this  author  and  his  associates  presented  passive  and  active  force  and 
moment  response  characteristics  of  major  human  joints,  associated  torques 
about  the  long-bone  axes  of  these  joints  and  some  aspects  of  joint 
modeling.  It  was  also  established  that  in  multisegmented  mathematical 
models  of  the  human  body  the  most  difficult  and  the  least  successful 
modeling  of  the  major  articulating  joint  has  been  the  shoulder  complex 
because  of  the  lack  of  appropriate  biomechanical  data  as  well  as  the 
anatomical  complexity  of  the  region. 

At  the  initiation  of  this  present  research  task,  a  partial  set  of 
data  existed  on  the  forced  kinematic  motion  (drawer  tests)  of  the  shoulder 
complex  of  one  male  subject.  Since  that  time  additional  data  were 
collected  on  two  more  male  subjects  to  provide  quantitative  results  on 
the  variability  of  the  "stiffness"  of  the  shoulder  complex  dependent  upon 
selected  orientations  of  the  upper  arm.  This  report,  first,  presents  a 
section  dealing  with  kinematics  by  means  of  sonic  emitters  and  a  special 
application  of  this  technique  to  the  shoulder  complex.  This  is  followed  by 
a  presentation  of  a  method  developed  to  analyze  the  raw  sonic  emitters 
data  by  utilization  of  the  three  most  accurate  sonic  emitters  out  of  six 
located  on  the  arm  cuff  of  the  test  subject.  The  final  numerical  results 
are  provided  in  the  form  of  plots  showing  the  passive  resistance  of  the 
shoulder  complex  versus  drawer  displacements  of  the  upper  arm  when  the 
humerus,  the  long-bone  of  the  upper  arm,  is  forced  axially. 


KINEMATICS  BY  MEANS  OF  SONIC  EMITTERS: 
APPLICATION  TO  THE  SHOULDER  COMPLEX 


The  quantitative  determination  of  the  nature  of  the  relative  motion 
between  two  body  segments,  which  are  connected  by  a  complex  anatomical 
joint,  is  of  prime  importance  to  the  biomechanician  as  well  as  to  those 
in  medicine.  The  simple  hinge  joint  with  one  degree  of  freedom  and  the 
ball  and  socket  joint  with  three  degrees  of  freedom  are  the  most  popular 
joints  employed  in  multisegmented  mathematical  models  of  the  human  body. 
Under  physiological  and  external  loads,  each  articulating  human  joint  can 
display  up  to  six  degrees  of  freedom  to  some  extent.  A  good  example  of  a 
general  joint  is  the  shoulder  complex  which  exhibits  four  independent 
articulations  among  the  humerus,  scapula,  clavicle,  and  thorax.  Of  course 
at  the  shoulder  complex,  the  six  degrees  of  freedom  refers  to  the  motion 
of  the  humerus  relative  to  a  fixed  portion  of  the  torso.  Since  flexing  of 
the  rib  cage  contributes  to  the  overall  range  of  motion  of  the  humerus, 
our  "fixed  body"  consists  of  the  lower  half  of  the  torso  (below  the  T-9 
vertebra)  and  the  upper  collateral  side  of  the  torso  (opposite  to  the 
shoulder  being  investigated).  The  kinematics  analysis,  which  is  considered 
here,  is  a  rigid  body  analysis.  It  is  very  essential  to  point  out,  before 
we  proceed  further,  that  no  human  body  segment  is  a  rigid  body  in  the  sense 
defined  in  mechanics.  However,  by  the  judicious  use  of  semi-rigid  orthotic 
shells  and  various  restraints,  approximate  rigid  body  conditions  were 
invoked  during  in-vivo  experiments  without  limiting  the  natural  range  of 
motion  of  the  joint  complex. 

In  this  study,  we  wish  to  monitor  the  displacement  of  a  point  in  the 
humeral  head  as  well  as  the  orientation  of  the  longitudinal  axis  of  the 
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humerus.  This  will,  in  essence,  give  us  some  idea  of  the  displacement  of 
the  glenohumeral  joint  with  respect  to  the  torso.  If  one  considers  the 
types  of  motion  possible  by  the  shoulder  complex  (i.e.  sweeping-type 
motions  as  well  as  drawer  type  motions)  it  becomes  obvious  that  the 
permissible  locations  of  the  glenohumeral  joint  define  a  volume  in  space 
relative  to  the  fixed  portion  of  the  torso.  In  addition,  the  voluntary 
boundaries  of  this  volume  are  extended  in  a  variable  manner  by  external 
loading  of  the  humerus.  The  role  of  the  sonic  digitizer  in  this  kinematics 
analysis  is  to:  1)  locate  the  fixed-body  origin  and  define  its  axis 
system  with  respect  to  the  microphone/sensor  assembly  coordinate  system, 
and  2)  to  monitor  the  motion  of  the  humerus  with  respect  to  that  same  axis 
system.  We  can  then  transform  the  humeral  kinematics  into  the  fixed-body 
system  of  the  torso  by  means  of  matrix  transformations. 

If  we  select  a  point  (origin)  in  the  moving  body  segment  and  establish 
a  local  axis  system  at  that  origin,  then  the  kinematics  of  the  moving  body 
segment  with  respect  to  the  fixed  reference  is  fully  defined  if  we  can 
continuously  monitor  the  location  of  the  selected  point  and  the  orientation 
of  the  local  axis  system  relative  to  the  fixed  reference  system.  This  can 
be  accomplished  by  tracking  three  points  (sonic  emitters)  fixed  on  the 
surface  of  the  moving  body  segment.  One  of  these  three  points  is  selected 
as  the  origin  and  unit  vectors  in  directions  from  the  origin  to  the  other 
two  points,  along  with  their  cross  product,  describe  the  local  axis  system. 

During  the  motion  analysis  it  is  desirable  to  know  the  positions  of 
the  three  emitters  simultaneously  at  any  given  instant.  Since  the  emitters 
are  fired  sequentially,  however,  this  was  not  possible.  One  way  to 


simulate  simultaneous  data  collection  would  be  to  use  a  "stop-frame"  type 
of  motion  in  which  the  moving  body  segment  would  be  motionless  during  any 
particular  "sweep"  of  the  three  emitters.  This  approach  was  ruled  out  for 
two  major  reasons:  1)  the  digitizer  cannot  take  commands  from  the 
computer,  any  pause  in  the  sequential  firing  order  makes  *cessary  that 
the  digitizer  is  reset  manually,  and  2)  the  necessity  of  ai  Ting  steadily 
increasing  loads  to  the  moving  body  segment  makes  this  tecT  ue 
undesirable.  A  more  tenable  solution  was  to  increase  the  l  .ig  rate  of 
the  digitizer  so  that  the  displacement  of  the  moving  body  segment  during 
any  given  sweep  of  the  three  emitters  would  be  small  compared  to  the 
distances  between  the  emitters  themselves. 

Two  major  modifications  of  the  data  collection  system  enabled  us  to 
maximize  the  digitizer  firing  rate.  First  an  IEEE-488  interface  in  place 
of  the  RS-232  interface  was  installed  between  the  GP6-3D  (sonic  digitizer) 
and  the  PDP-11/34.  While  the  RS-232  allows  single  line-serial  data 
transmission,  the  IEEE-488  uses  eight  lines  for  bit-parallel  byte-serial 
data  transmission.  Secondly,  the  programming  of  the  digitizer's  micro¬ 
processor  was  altered  so  that  slant  range  data  were  output  directly  to 
the  computer  (with  a  minimum  of  formatting)  thereby  eliminating  any  pre¬ 
processing  of  the  data  by  the  digitizer.  Although  these  changes  were 
very  time  consuming,  we  were  able  to  increase  the  emitter  pulse  rate  from 
15  pulses/sec  to  approximately  60  pulses/sec  (i.e.  &  20  sweeps  of  3 
emitters  each  second) . 

Initial  experiments  making  use  of  the  sonic  digitizing  technique 
demonstrated  that  additional  modifications  were  necessary  in  our  overall 


approach.  While  the  data  obtained  for  simple  planar  motions,  with  all  of 
the  emitters  facing  the  sensor  assembly,  are  quite  accurate;  more  general 
3-dimensional-type  motions  produced  data  of  more  questionable  accuracy. 

If  the  tip  of  an  emitter  is  rotated  away  from  the  sensor  assembly,  the 
sound  waves  must  travel  around  the  emitter1 s  base.  Likewise,  certain 
rigid  body  motions  either  partially  or  completely  block  the  emitter  from 
the  "view"  of  one  or  more  sensors.  These  conditions  can  either  aiter  the 
slant  range  data  or  produce  a  zero  reading  depending  on  the  amount  of 
blockage.  Therefore  it  was  necessary  to:  1)  devise  some  method  to  check 
the  accuracy  of  the  data  obtained,  and  2)  collect  redundant  data  in  order 
that  zero  readings  from  individual  emitters  would  not  affect  the  kinematics 
analysis.  The  following  section  addresses  these  considerations  and 
describes  a  new  data  collection  methodology  which  was  developed. 


A  NEW  DATA  COLLECTION  METHODOLOGY  5  DATA  ANALYSIS 

Since  the  kinematics  analysis  used  in  this  research  is  essentially 
based  on  vector  operations  involving  several  axis  systems,  the  following 
subscripts  are  utilized  to  specify  the  cartesian  coordinate  systems  and 
components  of  the  vectors  in  these  coordinate  systems : 

f:  fixed  body  system 

s:  sensor  assembly  axis  system 

l  or  (#,#,#):  arbitrary  or  specific  arm  cuff  local  axis  system, 
respectively. 


Likewise,  these  same  notations  are  used  as  superscripts  for  transformation 
matrices.  For  example  in  the  equation, 


As^  is  the  matrix  transformation  which  specifies  the  sensor  assembly  axis 

system  with  respect  to  the  fixed  body  system.  The  unit  vector  h  as  shown 

in  Fig.  2  defines  direction  of  the  long-bone  axis  of  the  humerus. 

During  the  tests,  the  subject  is  secured  in  a  chair  restraint  system 

which  restrains  the  lower  torso,  pelvis,  and  head  and  also  allows  the 

experimenter  to  position  the  subject  in  any  desired  orientation  relative 

to  the  sensor  assembly.  This  system  is  schematically  illustrated  in 

Fig.  1.  The  chair  has  motor-driven  pitch  and  yaw  capabilities  and  is 

adjustable  to  the  anthropometry  of  the  subject.  In  order  to  conveniently 

and  accurately  locate  the  fixed  body  system  relative  to  the  sensor 

assembly  axis  system,  a  relative  axis  locator  device  (RALD)  was  developed. 

The  RALD,  uses  four  sonic  emitters,  arranged  in  a  pyramidic  fashion,  to 

describe  an  axis  system  and  origin  location  in  space  relative  to  the  sensor 

assembly.  With  the  fixed  body  positioned,  as  desired,  with  respect  to  the 

sensor  assembly,  the  RALD  axes  are  aligned  with  the  fixed  body  system  and 

position  data  are  obtained  for  five  successive  sweeps  of  the  RALD’s  four 

emitters.  The  three  emitters  which  produce  the  most  consistent  readings 

—  f  s 

are  used  to  calculate  R  and  B  ,  where, 

s  —  * 

=  a  vector  describing  the  RALD  origin 
f  s 

EJ  =  a  transformation  matrix  defining  the  fixed  body  system  with 
respect  to  the  sensor  assembly  system 

Therefore : 
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Inversely,  the  sensor  assembly  axis  system  can  always  be  determined  with 
respect  to  the  fixed  body  system  by: 


There  is,  therefore,  no  need  to  measure  any  sensor  board  or  chair 

orientation  angles  as  was  done  previously.  This  approach  is  very  flexible 

from  a  data  acquisition  point  of  view.  Physical  dimensions  of  the  subject 

are  used  to  specify  a  vector  from  the  RALD  origin  to  the  actual  fixed  body 

—  —  s  f 

origin.  This  vector  F^.,  along  with  and  _B  are  stored  in  a  computer 

data  file  for  use  during  the  kinematic  analysis  program.  This  procedure  is 
easily  repeated  any  time  that  the  fixed  body  is  moved  with  respect  to  the 
sensor  assembly.  Next  we  address  the  application  of  sonic  emitters  to  the 
moving  body  segment,  i.e.  the  upper  arm. 

In  order  to  mount  emitters  on  the  upper  arm,  an  orthotic  cuff 
(schematically  shown  in  Fig.  2)  was  used.  This  cuff  consists  of  a  semi¬ 
rigid,  cylindrical,  plastic  shell,  which  is  encased  in  a  thin  rubber 
envelope.  It  extends  about  three-quarters  of  the  way  around  the  upper  arm. 
Velcro  straps,  at  the  top  and  bottom  of  the  cuff,  allow  it  to  be  tightly 
secured  circumferencially .  The  cuff’s  length  was  selected  based  on  the 
subjects’  upper  arm  dimensions  so  that  the  top  strap  rides  in  the  bicipital 
groove  (immediately  below  the  major  body  of  the  deltoid  muscle),  and  the 
bottom  strap  rides  in  the  antecubital  fold  (at  the  distal  end  of  the  biceps, 
above  the  elbow).  This  configuration,  along  with  the  static  friction  of 
the  rubber/skin  interface,  makes  the  cuff  suprisingly  secure  both 


3-5-6  AXIS 
SYSTEM 


LATERAL  VIEW  POSTERIOR  VIEW 


Fig.  2.  The  upper  arm  cuff,  its  sonic  emitters  and  one  of  the  local 
axis  systems  defined  by  emitters  3,5  and  6. 


longitudinally  and  circumferentially.  Six  sonic  emitters  are  affixed  to 
the  cuff  as  shown  in  Fig*  2.  Since  only  three  emitters  are  needed  to 
describe  the  location  and  orientation  of  the  upper  arm,  three  emitters  can 
be  totally  blocked  from  sensor  view  without  losing  kinematic  data.  A 
selection  of  the  ’’most  accurate"  three  emitters  is  made  in  cases  where 
four,  five,  or  all  six  emitters  produce  non-zero  data.  This  selection 
process  will  be  explained  shortly.  For  the  purposes  of  the  discussion  to 
follow,  we  will  denote  one  sweep  of  the  six  arm  cuff  emitters  by  the 
digitizer  as  a  record. 

Before  a  kinematic  test  is  performed,  an  initialization  procedure  is 
performed  to  establish  baseline  data  on  the  locations  of  the  six  cuff 
emitters  relative  to  the  humerus.  Let  us  consider  the  six  cuff  emitters 
as  points  on  the  upper  arm.  Referring  to  Fig.  2,  the  procedure  is  as 
follows.  The  elbow  is  flexed  to  90°,  and  one  emitter  is  positioned  above 
the  acromion  (pt.  7)  and  another  emitter  is  placed  below  the  olecranon 
process  of  the  ulna  (pt.  8).  A  data  file  which  contains  the  relative 
coordinate  locations  of  these  eight  emitters  is  collected.  A  unit  vector, 
h  ,  in  the  direction  from  pt.  7  to  pt.  8,  defines  the  long-bone  axis  of  the 
humerus.  Next,  knowing  the  coordinate  location  of  pt.  7,  we  use  hg  and  a 
pre-determined  length  based  on  the  subjects  dimensions,  to  calculate  the 
location  of  pt.  9.  The  point  9  is  an  approximation  of  the  center  of  a 
sphere  whose  surface  defines  the  articulating  surface  of  the  humeral  head. 
This,  therefore,  also  approximates  the  center  of  a  sphere  containing  the 
glenoid  cavity  surface  on  the  scapula.  The  point  9  is  what  we  define  as 
our  glenohumeral  joint  center,  whose  displacement  is  calculated  in  the 
kinematic  analysis . 


Turning  our  attention  to  the  six  arm  cuff  emitters,  we  note  that 
there  are  20  independent  combinations  of  three  emitters  (triads)  which  can 
be  used  to  describe  the  kinematics  of  the  upper  arm.  These  twenty  local 
systems  are  denoted  as  (a,b,c)  axis  systems,  where  a,  b,  and  c  refer  to 
the  cuff  emitter  numbers  which  are  used  in  the  grouping.  A  series- type 
representation  can  be  used  to  designate  the  numbered  sequences  for  the 
local  systems  by: 


(a,b,c) 


)) 


In  the  above  equation  the  summation-type  symbols  are  used  only  to  indicate 
how  the  notational  sequence  is  incremented,  and  are  not  meant  to  infer  any 
summation  operation  of  terms  as  in  a  series.  This  notational  scheme  insures 
that  all  six  arm  cuff  emitters  are  included  in  the  triad  groupings  and  that 
no  triad  i$  a  permutation  of  any  other  triad.  For  the  convenience  of 
ordering  with  respect  to  computer  calculations  the  triad  sequences  are 
labelled  1  through  20  according  to: 

Triad  #  Triad  Sequence 

1  (1,2,3) 

2  (1,2,4) 


20  (4,5,6) 


For  any  of  these  local  axis  systems,  emitter  b  is  the  origin,  a  unit  vector 
in  the  direction  from  a  to  b  describes  the  u-axis,  a  unit  vector  directed 


from  b  to  c  denotes  the  v-axis,  and  the  cross  product  of  u  x  v  is  w  which 
describes  the  w-axis.  Figure  2  shows  the  (3,5,6)  axis  system  as  an 
example.  Although  (3,5,6)  coordinate  system  is  an  orthogonal  system;  in 
general,  the  local  axis  systems  defined  by  the  triad  sequence  do  not 
necessarily  result  in  an  orthogonal  system.  For  example,  (1,3,6)  is  not 
an  orthogonal  coordinate  system  since  v  is  not  perpendicular  to  u. 

A  A  a 

However,  it  is  an  easy  task  to  define  a  new  v  such  that  v^Lto  u  by 
shifting  the  origin  of  the  coordinate  system  along  the  unit  vector  u.  Th 
all  local  axis  systems  become  an  orthogonal  one.  Note  that  the  x,  y,  and 
z  components  of  u,  v,  and  w  compose  a  direction  cosine  matrix  given  by: 


l(3,5,6)s 


v 

z 


w 

z-l 


This  matrix  defines  the  local  (3,5,6)  system  relative  to  the  x,y,z  axis 
system  of  the  sensor  assembly.  Referring  again  to  Fig.  2  if  we  now 
calculate  the  vector  from  pt.  5  to  pt.  9,  (Jg) ,  we  can  then  locate  the 
joint  center  with  respect  to  the  local  system  by: 


J(3,5,6) 


l(3,5,6)s  j 


and  likewise: 


h(3,5,6) 


l(3,5,6)s  ~ 


In  the  initialization  computer  program,  vectors  locating  the  joint 
center  and  describing  the  humeral  long  axis  are  calculated  for  all  twenty 


axis  systems  and  this  information  is  placed  in  a  data  file.  The 
initialization  program,  however,  also  establishes  our  "accuracy  criteria" 
for  local  triad  selection.  One  must  realize  that  during  a  kinematic  test 
we  do  not  know  ahead  of  time  which  three  emitters  will  provide  the  most 
accurate  description  of  the  rigid  body  location/orientation  for  any  given 
record.  During  the  initialization  procedure,  however,  we  establish  a 
"best  condition"  data  set  that  records  the  relative  positions  of  the  six 
cuff  emitters.  We  then,  in  turn,  use  these  data  to  calculate  the 
relative  orientations  and  locations  of  the  twenty  local  axis  systems  with 
the  arm  cuff/arm  in  optimal  view  of  the  sensor  assembly. 

For  any  two  axis  systems  (i.e.  the  i-th  and  j-th  axis  systems,  where 
i  and  j  are  between  1  and  20)  we  calculate  the  transformation  matrices 
with  respect  to  the  sensor  assembly  axis  system: 


Therefore,  the  relationship  of  the  j-th  system  with  respect  to  the  i-th 
system  is  given  by 


where  T1J  is  3x3  transformation  matrix  describing  the  i-th  system  with 
respect  to  the  j-th  system.  The  transformation  matrices  relating  each  of 
the  20  local  systems  to  every  other  system  are  calculated  in  this  manner, 
thereby  arriving  at  190  direction  cosine  matrices.  Only  transformations 
in  one  direction  are  calculated  since  it  is  obvious  that: 


\ 

J 

J 

J 

i 


« 

i 


i 

i 
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These  matrix  relationships  are  stored  in  the  same  data  file,  along  with 
the  joint  center  vectors  and  the  humeral  axis  vectors  for  the  twenty 
local  systems.  The  final  step  in  the  initialization  procedure/program  is 
to  calculate  the  distances  between  the  various  local  axis  systems  and  to 
store  these  data  along  with  the  previously  mentioned  information.  Note  that 

these  initial  or  baseline  data  are  obtained  with  the  arm  cuff  emitters  and 
humeral  axis  emitters  in  optimal  view  of  the  sensor  assembly. 

For  a  typical  test,  with  the  upper  arm  in  motion,  the  choice  of  the 
,fmost  accurate”  local  system  is  made  during  the  kinematic  analysis  on  a 
record-by-record  basis.  For  each  record  of  kinematic  data,  the  coordinate 
data  of  the  six  cuff  emitters  are  once  again  used  to  calculate  the  direc¬ 
tion  cosine  matrices  describing  the  local  axis  systems  with  respect  to  the 
sensor  assembly  system.  These  new  (kinematic)  local  axis  systems  are,  in 
turn,  used  to  obtain  the  matrix  transformations  (interrelationships) 
between  the  local  axis  systems.  Let  us  designate  these  matrix 

transformations  with  T^.  ^  .  . 

—  kinematic 

The  transformation  matrices  obtained  from  the  kinematic  data  will 
differ  from  their  initialized  counterparts  for  the  following  reasons: 

1)  Since  the  emitters  do  not  fire  simultaneously,  the  "apparent" 
axis  system  for  a  given  emitter  triad  utilizes  data  obtained  at 
three  different  times. 

2)  As  emitters  are  rotated  away  from  the  "view"  of  the  sensor 
assembly,  the  acoustic  pulse  may  still  be  detected,  but  the 


distance  measurements  are  dilated  because  the  acoustic  wave 


needs  to  travel  around  part  of  the  emitter  itself. 

3)  Acoustic  blockage  similar  to  item  2  may  be  caused  by  the 
location  of  part  of  an  object  such  as  a  body  segment  between  the 
emitters  and  the  sensor  assembly.  Obviously,  the  object  affects 
the  acoustic  field  which  is  produced  by  the  emitter  pulse  and 
will  cause  dilations  in  the  distance  measurements. 

4)  The  reproducibility  of  the  measurements  varies  ±  1  in  the  last 
significant  digit  which  corresponds  to  ±  .1  mm  or  ±  • 1%  of  the 
distance  measured. 

Therefore,  during  data  analysis,  it  is  the  deviation  from  the 
predetermined,  initialized  axis  system  interrelationships  that  constitutes 
the  selection  procedure  for  establishing  the  local  axis  system  which  best 
describes  the  moving  body  segment  location  and  orientation.  If  there  were 
no  errors  in  measurements,  and  a  ’’stop-motion”  type  of  sweeps  of  the  cuff 
emitters  were  possible,  then  we  should  obtain  the  equality: 


kinematic  initial 


or 


(iij  )(iij  )_1=[i] 

\  kinematic  /  \  initial  / 

where  [I]  is  the  3x3  identity  matrix.  This,  however,  is  not  the  case,  and 
a  general  matrix  [G]  with  off-diagonal  terms  is  obtained.  Therefore: 
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The  matrix  [G] ,  thus,  represents  an  "apparent”  rotation  of  the  i-th  axis 
system  with  respect  to  the  j-th  axis  system.  The  angle  of  rotation 
between  these  two  coordinate  systems  is  given  by  [24]: 


Y-  •  =  cos 
ij 


‘(1C 


gll  +  g22  +  g33 


■1]) 


This  "apparent11  or  perceived  rotation,  which  results  from  the  measurement 
error  factors  previously  mentioned  is  denoted  as  the  skew  angle  error, 
y^y  Note  that  the  skew  angle  error  will  affect  joint  center  calculations. 
If  we  use  an  average  distance  of  20  cm  from  a  local  axis  system  to  the 
joint  center  and  a  skew  angle  error  of  0.20  degrees,  the  magnitude  of  the 
maximum  error  of  in  the  joint  center  computation  will  be 

(6 i j ) sk  =  20  sin (0.20)  =  0.070  cm 

In  this  way,  (6  —  )^  values  are  calculated  for  all  190  matrix  inter¬ 
relationships.  If  we  likewise  compare  the  distances  between  origins  for 
th 3  kinematic  data  record  versus  the  initialized  data  we  can  obtain  an 
error  that  can  be  induced  in  the  joint  center  calculation  due  to  apparent 
local  axis  system  dislocation,  (6.^)^,  where: 

(6  i  .)d  =  D..  -  IK. 

1Jkinematic  initial 


with  D. . 


and  D. . 


^ initial 


as  the  intra-origin  distance  between  the 


^kinematic 

i-th  and  j-th  system  origins  for  the  kinematic  and  initialized  records, 
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respectively.  Having  now  calculated  (6  —  )^  and  (6  —  )^  f°r  the 

local  axis  system  interrelationships,  a  least-squares  analysis  is 
performed  to  select  the  local  axis  system  which  minimizes  the  interrela- 
tional  errors.  For  each  i-th  system  a  (<5-^)avg  is  obtained  from: 


20 

OOavg  '  15  l(6ii)sk  *  N  ) 


2( 

d  \ 


1/2 


for  i  /  j 


From  a  probability  standpoint,  the  system  having  the  smallest  C^)avg  is 

likely  to  be  the  most  accurate  with  respect  to  the  sensor  assembly.  A 

physical  interpretation  is  that  the  system,  thus  selected,  perceives  the 

least  amount  of  change  among  the  systems  surrounding  it;  likewise  systems 

that  perceive  large  changes  in  many  surrounding  systems  are  more  likely  to 

be  in  error  themselves.  Having  thus  selected  the  "best-fit"  triad  of 

emitters  from  the  six  arm  cuff  emitters  we  denote  the  axis  system  that  they 

describe  with  respect  to  the  sensor  assembly  by  means  of  the  transformation 
.Sts  J£s 

matrix  L_  .  The  transformation  JL  is  used,  along  with  the  retrieved 
initialized  data  on  the  joint  center  vector  and  humeral  axis  vector  for  that 
particular  axis  system,  to  help  calculate  the  joint  center  location  and 
humeral  axis  orientation  with  respect  to  the  fixed  body  segment. 

Referring  to  Fig.  3,  the  vector  to  the  origin  of  the  selected  local 
system  on  the  arm  cuff  is  denoted  as  C^.  By  inspection,  the  vector  from 
the  fixed  body  origin  to  the  joint  center  (Gf)  can  be  calculated  by: 

(G)f  =  (-F-R  +  C  +  J)f 


where  F^.  is  given  in  the  fixed  body  system,  and: 


r*VV7V 


Cc  =  Bfs  C 

f  —  S 

Jf  =  Bfs(L£s)T  Jz 


This  procedure  is  repeated  for  each  record;  each  time  a  new  joint  center 
location  and  humeral  axis  orientation  are  computed  with  respect  to  the  fixed 
body  segment  (the  torso).  Since  these  records  determine  the  instantaneous 
positions  and  orientations  of  the  moving  body  segment,  the  upper  arm,  the 
kinematics  of  the  moving  body  segment  is  thus  fully  known. 

To  complete  the  data  analysis  one  must  describe  the  function  of  the 
force  applicator  which  can  be  seen  in  Fig.  4  showing  a  subject  in  position 
for  superior- inferior  drawer  testing.  The  force  applicator  has  two 
functions:  (1)  to  determine  the  direction  and  the  magnitude  of  the  force 

application,  and  (2)  to  determine  the  location  of  the  force  application  on 
the  moving  body  segment  with  respect  to  the  sensor  assembly  axis  system. 

Th  is  is  accomplished  by  monitoring  three  strategically  located  sonic  emitters 
on  the  force  applicator.  The  force  applicator  actually  has  four  emitters  (no.  1 
is  located  on  the  opposite  side  from  no.  3),  but  only  three  are  in  view 
of  the  sensor  assembly  for  any  given  test.  For  the  forced  motion  tests, 
the  subject  and  sensor  assembly  are  arranged  in  such  a  way  as  to  provide 
optimum  tracking  of  the  force  applicator.  For  this  reason,  redundant 
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emitters,  as  in  the  case  of  the  arm  cuff,  were  not  necessary.  Since  the 
exact  relative  positions  of  the  three  sonic  emitters  are  known,  sufficient 
information  exists  for  a  vector  analysis  to  determine  both  the  location 
and  direction  of  the  applied  force. 

If  we  refer  to  Fig.  5,  we  see  that  the  arrangement  of  the  emitters 
allows  us  to  readily  specify  the  orthogonal  axis  system  of  the  transducer. 
Initially,  a  unit  vector  in  a  direction  from  pt.  1  to  pt.  2  defines  the 
z-axis  of  the  transducer.  Next,  pt.  4  is  calculated  as  the  mid-point 
between  points  2  and  3.  Let  dj  specify  the  magnitude  of  the  vector  from 
pt.  2  to  pt.  3,  then  the  vector  is  obtained  by: 

L,  =  (  T  (U12  X  U23>  )s 

where  the  subscript  s,  as  before,  denotes  that  the  vectors  in  the  above 
equation  are  with  respect  to  the  sensor  assembly  axis  system.  Using  this 
same  convention,  we  denote  the  position  vector  which  locates  the  number 
2  emitter  as  (.r2^sm  Therefore,  by  inspection  of  Fig.  5,  we  see  that  the 
location  of  pt.  5  (on  the  longitudinal  axis  of  the  force  applicator)  is 
given  by: 


(r2>s  +  T  <-u23>s  + 


P 

s 


and  the  location  of  pt.  6  (the  point  of  force  application)  is  given  by: 


Note  that  unit  vectors  from  pt.  5  to  points  2  and  3  define  the  transducer's 
x  and  y  axes,  respectively,  with  respect  to  the  sensor  assembly  axis  system 
Since  both  the  point  of  force  application  and  the  transducer  axes  are  now 
known  with  respect  to  the  sensor  assembly  system,  these  data  are  easily 
transformed  into  the  fixed  body  system  as  was  done  with  the  arm  cuff  data. 
Once  again,  the  data  are  updated  on  a  record  by  record  basis. 

RESULTS 

Using  the  concepts  presented  in  the  previous  section  the  data  obtained 
from  a  number  of  tests  which  were  performed  for  determination  of  the 
passive  resistance  of  the  shoulder  complex  to  forced  (drawer- type)  motions 
were  analyzed.  In  these  tests  the  upper  arm  was  forced  axially  for  a 
variety  of  upper  arm  orientations  with  respect  to  torso.  Before  we  present 
the  results  of  the  data  analysis  we  will  briefly  describe  the  testing 
procedure . 

During  the  tests,  the  subject  is  secured  in  the  chair  restraint 
system  which  allows  unimpeded  movement  of  the  shoulder  complex.  After 
the  initialization  procedure  is  completed,  and  the  fixed  body  is  located 
with  respect  to  the  sensor  assembly  by  means  of  the  RALD,  the  force 
applicator  is  fixed  to  the  arm  by  means  of  a  rigid  cuff  extension  so  that 
the  force  can  be  applied  comfortably  along  the  humeral  axis.  In  addition, 
the  rigid  cuff  places  the  upper  arm  in  a  position  of  approximate ly  25c  of 
medial  rotation.  This  position  is  near  the  mid-range  of  the  shoulder 
medial- lateral  rotation  limits  and  is  a  comfortable  position  for  the 
subject. 


For  each  upper  arm  orientation  defined  in  the  next  paragraph,  the 
arm  is  forced  proximally  and  then  distal ly  along  its  longitudinal  axis. 
Before  each  test,  the  subject  was  asked  to  allow  his  shoulder  complex  to 
assume  its  most  comfortable,  "natural"  position.  This  established  the 
neutral  position  for  each  drawer.  During  the  tests  the  subjects  were 
instructed  to  avoid  muscular  resistance,  and  each  test  was  terminated  when 
the  subject  expressed  discomfort.  Three  male  subjects  were  tested*  at  six 
positions  of  the  upper  arm  orientations.  Subjects  were  university  students 
with  no  special  training  in  athletics.  Selected  anthropometric  measure¬ 
ments  of  the  subjects  are  given  in  Table  I.  The  definitions  of  the 
anthropometric  measurements  provided  in  Table  I  can  be  found  in  [3  §  13]. 
The  six  positions  can  be  defined  by  means  of  0  and  (j>  angles.  The  6  angle 
refers  to  the  angle  between  the  z-axis  of  the  torso  and  the  long-bone  axis 
of  the  upper  arm;  this  angle  also  defines  the  shoulder  flexion-extension 
in  the  sagittal  plane.  The  <j>  angle  refers  to  the  angle  between  the 
projection  of  the  long-bone  axis  of  the  upper  arm  on  the  xy-plane  and  the 
x-axis;  the  positive  and  negative  values  of  this  angle  also  define  the 
shoulder  abduction  and  adduction,  respectively. 

The  test  data  were  analyzed  and  presented  for  the  following  six 
positions  of  the  upper  arm:  1)  0  =  0°,  <f>  =  0°  which  corresponds  to  the 
inferior- superior  drawer  test  as  shown  in  Fig.  4;  2)  0  =  90°,  <J>  =  0° 
which  corresponds  to  anterior-posterior  drawer  test  in  the  sagittal  plane; 
3)  0  =  90°,  <p  =  30°;  4)  6  =  90°,  <J>  =  60°;  5)  0  =  90°,  <j>  =  90°  which 

*  The  testing  was  conducted  in  conformance  with  the  Human  Subject  Program 
Guidelines  established  by  the  Ohio  State  University  Human  Subject 
Review  Comnittee. 
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TABLE  I 


SELECTED  ANTHROPOMETRY  OF  SUBJECTS  1,  2,  3 
(ALL  LENGTH  DIMENSIONS  ARE  IN  CM) 


12  3 


1 

WEIGHT  (NEWTONS) 

787.6 

851.9 

816.1 

STATURE 

188.6 

185.5 

186 

SHOULDER  CIRCUMFERENCE 

115 

126.4 

116.8 

WAIST  CIRCUMFERENCE 

96 

101.6 

99 

WRIST  CIRCUMFERENCE 

17 

17.8 

17.6 

LOWER  ARM  CIRCUMFERENCE 

29.8 

29.8 

28.9 

BICEPS  CIRCUMFERENCE 

32.2 

37.1 

32 

THIGH,  UPPER  CIRCUMFERENCE 

57.5 

57.1 

57.5 

THIGH,  LOWER  CIRCUMFERENCE 

43 

40.6 

44.5 

CALF  CIRCUMFERENCE 

41 

40.6 

38 

ANKLE  CIRCUMFERENCE 

26 

27 

26.5 

FOREARM-WRIST  LENGTH 

23.5 

21.5 

21.5 

SHOULDER-ELBOW  LENGTH 

41 

38.5 

38.5 

SHOULDER  HEIGHT,  SITTING 

67.5 

67 

68.5 

SITTING  HEIGHT 

97.2 

97 

94 

SHOULDER  BREADTH _ 50 _ 50.5 _ £7 

CHEST  BREADTH  34.5  35  33 


CHEST  DEPTH 

23 

23 

22 

WAIST  DEPTH 

20.5 

22.5 

21 

BUTTOCK- KNEE  LENGTH 

66 

62 

64  i 

| 

BUTTOCK-POPLITEAL  LENGTH 

54.5 

50.8 

53.5  ! 

KNEE  HEIGHT,  SITTING 

61 

58.5 

61.5 

ELBOW-TO-ELBOW  BREADTH 

46.5 

48 

47.5 

HIP  BREADTH  SITTING _ 41 _ 37 

KNEE-TO-KNEE  BREADTH,  SITTING  25  25.5 


38 _ 

22.5 


V  -S 


corresponds  to  medio- lateral  drawer  test  in  the  frontal  plane;  and 
6)  9  =  90°,  <(>  =  115°.  The  numerical  results  for  these  are  presented  in 
the  form  of  plots  in  Figs.  6-23  showing  the  passive  resistance  of  the 
shoulder  complex  as  a  function  of  drawer  displacement  of  the  upper  arm 
along  its  long-bone  axis.  In  these  plots  the  ordinate  represents,  in 
Newtons,  the  magnitude  of  the  axially  applied  force  and  the  abscissa 
represents,  in  centimeters,  the  computed  components  of  the  joint  center 
displacement  along  the  force  application  direction.  The  force  and 
displacement  values  were  computed  on  a  record  by  record  basis.  In 
addition,  for  each  record,  the  average  skew  error  and  average  intra-origin 
error  for  the  local  axis  system  that  was  selected  to  describe  the  upper 
arm  position  were  computed,  and  displayed  on  the  computer  printout.  These 
average  errors  pertain  to  the  "apparent"  errors  in  all  the  i-th  local 
systems  with  respect  to  the  j-th  system  which  is  selected  for  use  in  joint 
center  calculations.  For  the  presented  results,  the  average  skew  error 
was  less  than  .45  degrees  and  the  average  intra-origin  error  was  less  than 
.15  centimeters. 
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MAGNITUDE  OF  RESISTIVE  F 


NEUTRAL  POSITION  (cm) 


Fig,  6.  Passive  resistive  force  versus  drawer  displacement  for  the 
superior-inferior  drawer  test  of  the  shoulder  complex  of  the  first 
subject  for  the  upper  arm  orientation  along  torso  (9  =  O'  ,  <J>  =  0°)  . 


y  300 

or 

o 

Ll 

250 


h 

tI  200 

CO 

LlI 

or 


100 


1 

t 

— 

t  SUPERIOR 

INFERIOR _ . 

:  DRAWER 

DRAWER  j 

i  i  i° 

i  i  i 

10  864202468  10  12 

DRAWER  DISPLACEMENT  FROM 
NEUTRAL  POSITION  (cm) 


Fig,  7.  Passive  resistive  force  versus  drawer  displacement  for  the 
superior-inferior  drawer  test  of  the  shoulder  complex  of  the  second 
subject  for  the  upper  arm  orientation  along  torso  (0  =  0°,  <J>  =  0°). 
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Fig.  8.  Passive  resistive  force  versus  drawer  displacement  for  the 
superior- inferior  drawer  test  of  the  shoulder  complex  of  the  third 
subject  for  the  upper  arm  orientation  along  torso  (0  =  0°,  <J>  =  0°). 
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Fig.  9.  Passive  resistive  force  versus  drawer  displacement  for  the 
anterior-posterior  drawer  test  of  the  shoulder  complex  of  the  first 
subject  (0  =  90°,  <J>  =  0°). 
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Fig.  11.  Passive  resistive  force  versus  Jr.  wer  displacement  for  the 
anterior-posterior  drawer  test  of  the  shculu^r  complex  of  the  third 
subject  (9  =  90°,  <j>  =  0°). 
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Fig.  12.  Passive  resistive  Force  versus  drawer  displacement  of  the 
shoulder  complex  of  the  first  subject  for  the  0  =  90°,  <f>  =  30°  upper 
arm  orientation. 
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Fig.  13.  Passive  resistive  force  versus  drawer  displacement  of  the 
shoulder  complex  of  the  second  subject  for  the  0  =  90c,  =  50°  upper 

arm  orientation. 
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Fig.  14.  Passive  resistive  force  versus  drawer  displacement  of  the 
shoulder  complex  of  the  third  subject  for  the  0  =  90°,  <f>  =  30°  upper 
arm  orientation. 
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Fig.  15.  Passive  resistive  force  versus  drawer  displacement  of  the 
shoulder  complex  of  the  first  subject  for  *  •  0  =  90°,  <j>  ~  60°  upper 

arm  orientation. 
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Fig.  16.  Passive  resistive  force  versus  drawer  displacement  of  the 
shoulder  complex  of  the  second  subject  for  the  0  =  90°,  <j>  =  60°  upper 
arm  orientation. 
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Fig.  17.  Passive  resistive  force  versus  drawer  displacement  of  the 
shoulder  complex  of  the  third  subject  for  the  8  =  90°,  <f>  =  60°  upper 
arm  orientation. 


MAGNITUDE  OF  RESISTIVE  FO 


MAGNITUDE  OF  RESISTIVE  FORCE  (N) 


350 

300 

250 

200 

150 

100 

50 

0 

10  8  6  4  2  0  2  4  6  8  10  12 

DRAWER  DISPLACEMENT  FROM 
NEUTRAL  POSITION  (cm) 


Fig.  20.  Passive  resistive  force  versus  drawer  displacement  of  the 
shoulder  complex  of  the  third  subject  for  the  0  =  90°,  <f>  =  90°  upper 
arm  orientation. 
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Fig.  21.  Passive  resistive  force  versus  drawer  displacement  of  the 
shoulder  complex  of  the  first  subject  for  the  0  =  90°,  $  =  115°  upper 
arm  orientation. 
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Fig.  22.  Passive  resistive  force  versus  drawer  displacement  of  the 
shoulder  complex  of  the  second  subject  for  the  9  =  90°,  <J>  =  115°  upper 
arm  orientation. 
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Fig,  23.  Passive  resistive  force  versus  drawer  displacement  of  the 
shoulder  complex  of  the  third  subject  for  the  0  =  90°,  <J>  =  115°  upper 
arm  orientation. 


CONCLUDING  REMARKS 


The  plots  of  the  forced  kinematic  motion  (drawer  tests)  of  the 
shoulder  complex  presented  in  this  report  provide  some  quantitative 
results  on  the  variability  of  the  "axial  stiffness"  of  the  human  shoulder 
complex  dependent  upon  the  upper  arm  orientation.  Additional  quantitative 
assessment  of  the  axial  stiffness  of  the  shoulder  complex  can  be  made  by 
curve  fitting  the  numerical  results  plotted  in  Figs.  6-23.  Two  functions 
in  the  form  of  y^  =  ae^X  and  y ^  =  cx^  are  chosen.  Considering  the  fact 
that  both  exponential  and  power  forms  are  the  two  functional 
representations  which  predominate  in  the  elastic  models  for  the  constitutive 
stress-strain  behavior  of  fibrous  connective  tissues,  choice  of  these  two 
functions  is  quite  natural.  The  exponential  form  was  originally  proposed 
by  Fung  [19]  and  later  used  in  various  versions  [20,21,28,29,30];  the 
usage  of  power  form  was  equally  wide-spread  [22,23,25,26,27].  It  is  also 
important  to  point  out  that  the  axial  stiffness  characteristics  of  the 
shoulder  complex  are  closely  related  to  the  constitutive  behavior  of  the 
soft  tissues  associated  with  the  shoulder  complex. 

Table  II  contains  the  coefficients,  a,  and,  b,  for  the  exponential 

curve.  The  numbers  in  parentheses  in  both  Tables  II  and  III  designate  the 

2 

regression  coefficients  of  determination,  r  ,  which  indicate  the  degree  of 
fit  attained  with  the  tabulated  coefficients.  Variation  of  the  b 
coefficient  for  the  exponential  function  with  respect  to  upper  arm 
orientation  is  shown  in  Fig.  24  for  two  subjects.  The  coefficients,  c, 
and,  d,  for  the  power  curve  fitting  are  given  in  Table  III.  In  certain 
cases  the  power  curve  appears  to  be  the  more  appropriate  choice.  In  the 


curve  fitting  process  for  both  types  of  functions  the  data  points  near 
the  neutral  position  are  excluded.  In  this  region  the  curves  exhibit  a 
rapid  increase  initially,  followed  by  a  region  of  less  rapid  increase.  It 
is  believed  that  this  phenomenon  is  an  artifact  of  the  testing  procedure. 
It  is  possibly  due  to  a  reflex  resistance  caused  by  the  activation  of 
stretch  receptors  in  the  muscles  of  the  shoulder  complex. 

Finally,  although  there  are  some  intra-subject  variations  in  the 
behavior  of  the  axial  stiffness  of  the  shoulder  complex  for  the  subjects 
tested,  strong  similarities  in  the  results  do  exist  and,  thus,  certain 
characteristic  response  patterns  for  the  axial  stiffness  of  the  human 
shoulder  complex  can  be  identified. 
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TAB  LI:  II 


COi;m Cl  l:\i  s.  a  and  b,  OF  THE  EXPONENTIAL  CURVE,  y  =aeL 


UPPER  ARM  ORIENTATION 
ti  TYPE  OF  DRAWER  TEST 


SUPERIOR  DRAWER 


INFERIOR  DRAWER 


ANTERIOR  DRAWER 


=  90  POSTERIOR  DRAWER 


30  OUTWARD  DRAWER 


=  90  INWARD  DRAWER 


=  60  OUTWARD  DRAWER 


90w  INWARD  DRAWER 


=  90  LATERAL  DRAWER 


90  MEDIAL  DRAWER 


=  115  OUTWARD  DRAWER 


=  90  INWARD  DRAWER 


1 - 

SUBJECT 

1 

a 

b 

12 

602 

0. 

279 

10 

999) 

8 

248 

0. 

531 

(0 

999) 

31 

457 

0. 

102 

(0 

987) 

15 

040 

0. 

183 

(0 

991) 

4, 

.008 

0. 

259 

(0 

.999) 

7 , 

.209 

0. 

173 

(0, 

.999) 

16. 

.855 

0. 

241 

(0. 

.998) 

15. 

.614 

0. 

178 

(0. 

.997) 

20. 

.  159 

0. 

521 

(0 

.994) 

54.515  0.185 

10.907 j 


SUBJECT  2 
a  b 

7.191  0.569 

(0 .999) 

6.538  0.553 

(0.992) 

4.619  0.207 

(0.995) 

8.672  0.209 

(0.993) 

3.763  0.264 

(0.995) 

1.006  0.319 

(0.971) 

3.602  0.227 

(0.969) 

1.253  0.364 

(0.998) 


5.075  0.305 

(0.978) 


26.588  0.54’  :  12.092  0.410 

(0.988:  j  (0.996) 

29.965  0.540  [  15.508  0.452 

(0.993)  j  (0.988) 


SUBJECT  3 
a  b 

20.171  0.224 

(0.972) 

11.179  0.280 

(0.996) 

12.656  0.172 

(0.998) 

4.656  0.210 

(0.997) 

20.754  0.152 

(0.994) 

9.989  0.166 

(0.998) 

16.386  0.403 

(0.999) 

19.469  0.153 

(0.972) 

27.593  0.365 

(0.993) 

39.015  0.150 

(0.999) 

40.437  0,285 

(0.985) 

16.560  0.328 

(0.999) 


TAB  LI:  111 


COEFFICIENTS,  c  and  d,  01  Till:  POWER  CURVL,  y,  =cx 


UPPER  ARM  OR  I  ENTAT  I  ON 
I,  TYPF;  OF  DRAW  HR  ThST 


=  0  SUPERIOR  DRAWER 


=  0  INFERIOR  DRAWER 


4>  =  0  ANTERIOR  DRAWER 


0  =  90  POSTERIOR  DRAWER 


30  OUTWARD  DRAWER 


SUBJECT  1 


SUBJECT  2 


1.707  2.071  I  2.230  1.939 

(0.992)  (0.991) 

0.554  2.617  1.092  2.268 

(0.983)  (0.971) 

12.015  0.949  0.137  2.470 

(0.920)  (0.961) 

1.437  1.849  1.340  1.766 

(0.999)  (0.970) 

0.120  2.699  0.118  2.698 


(0.984) 


(0.967) 


0.359  2.102  0.003  3.866 

(0.983)  (0.948) 

12.920  1.083  0.090  2.633 

(0.939)  (0.925) 


SUBJECT  3 


4 . 334  1 . 659 

(0.999) 

1.505  2.121 

(0.980) 

3.002  1.477 

(0.941) 

0.105  2.601 

(0.975) 

6.407  1.231 

(0.967) 


18.723  1. 

(0.952) 


LATERAL  DRAWER 


MEDIAL  DRAWER 


OUTWARD  DRAWER 


INWARD  DRAWER 


7.339  1.709 

(0.996) 

16.080  1.071 

(0.928) 

20.780  1.244 

(0.998) 

29.043  1.136 

(0.986) 


0.551  2.309 

(0.936) 

3.325  2.126 

(0.991) 

8.121  1.825 

(0.999) 


15.461  1. 

(0.977) 

27.260  1. 

(0.991) 

5 . 39  7  1  . 

(0.984) 


COEFFICIENT, 


Rg.  -4.  The  coefficients,  b,  of  the  exponential  curve  fitting  of  the 
passive  resistance  data  of  the  shoulder  complex  at  various  upper  arm 
positions  (For  the  location  indicated  by  ?  data  are  not  available). 
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